Quantum frequency estimation with trapped ions and atoms 
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We discuss strategies for quantum enhanced estimation of atomic transition frequencies with ions 
stored in Paul traps or neutral atoms trapped in optical lattices. We show that only marginal 
quantum improvements can be achieved using standard Ramsey interferometry in the presence of 
collective dephasing, which is the major source of noise in relevant experimental setups. We therefore 
analyze methods based on decoherence free subspaces and prove that quantum enhancement can 
readily be achieved even in the case of significantly imperfect state preparation and faulty detections. 
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The ultra-precise estimation of physical parameters is 
of great importance for countless applications, such as 
atomic clocks, gravitational wave detectors, laser gyro- 
scopes or microscopy. Quantum enhanced precision mea- 
surements have the potential to significantly increase the 
precision of parameter estimation compared to classi- 
cal methods in such applications This is generally 
achieved by preparing a quantum probe-state which has 
a higher sensitivity with respect to the quantity to be 
estimated. Assuming that this probe consists of non- 
interacting, identical subsystems (e.g. N particles), the 
estimation uncertainty in many applications including 
optical or atomic interferometry can then ideally be im- 
proved from the standard quantum limit (SQL), which 
scales like ^h/N, to the Heisenberg limit, which scales 
like 1/^ @,|3l. Endeavors to attain the Heisenberg limit 
(or at least to beat the SQL) are made in many branches 
of physics including quantum photonics Q and atomic 
physics [sl-fioj. In atomic physics, the measurement of 
atomic transition frequencies with Ramsey interferome- 
try has been established as an important tool, not only 
for general spectroscopic purposes but also to determine 
frequency standards on which atomic clocks are based 
on Improvements of Ramsey interferometry via 

quantum effects are therefore highly desirable. As in 
other quantum technologies like quantum computing and 
communication, the biggest obstacle for the realization 
of such a quantum interferometer is the presence of un- 
avoidable noise and imperfections. A practical quantum 
sensor must therefore use probe-states which are robust 
under realistic circumstances, as well as preparation and 
detection schemes which can be performed with high fi- 
delity. 

In this paper we analyze methods for quantum 
enhanced estimation of atomic transition frequencies 
with Ramsey interferometry, and generalizations thereof, 
which can improve the measurement uncertainty to the 
Heisenberg limit in the presence of noise, and which toler- 
ate imperfect state preparation and detection. A scheme 
for quantum enhanced Ramsey interferometry has been 
proposed some time ago [l^, but it has subsequently 
been shown that in the presence of noise, in form of 



uncorrelated dephasing, the scheme has only little or 
no advantage compared to its classical counterpart [l3j . 
However, recent experimental breakthroughs with closely 
spaced particles, particularly ions stored in linear Paul 
traps [l4l - [l7f . show that the major source of noise in 
these systems consists of correlated dephasing. Moti- 
vated by this insight we first analyze conventional Ram- 
sey interferometry and show that in the presence of cor- 
related dephasing hardly any quantum enhancement can 
be achieved. We therefore discuss alternative methods 
which make use of decoherence free subspaces [l^ and 
show that they lead to quantum enhanced precision even 
in the presence of significantly imperfect state prepara- 
tion and faulty detections. Our approach is mainly mo- 
tivated by recent experiments with trapped ions, but it 
can also be applied to cold atoms stored in optical lat- 
tices [l^. The main body of this paper concisely sum- 
marizes our results. Details of calculations can be found 
in the appendices. 

We consider N two-level atoms or ions, with internal 
states |0) and |1), which are, e.g., stored in a linear Paul 
trap. These atoms are prepared in an input state \ipin) 
and undergo the process shown in Fig. [IJa), i.e. each 
atom accumulates a phase ip during a time t and is fi- 
nally measured. The process is repeated v times and 
based on the measurement outcomes the phase can be 
estimated. For simplicity we assume that the measure- 
ment and preparation times are much smaller than t, such 
that the total time of the experiment is given by T = i^t. 
It is our goal to make the uncertainty of the estimated 
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FIG. 1: (a) Schematic Ramsey interferometer, (b) General- 
ized setup where atoms acquire different phases (see text). 
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phase as small as possible for a given T and N . 

In conventional Ramsey interferometry the input state 
is given by the product state l^-fD = [(|0) + \l))/V2]®^ 
which is prepared by a 7r/2-pulse using a laser with fre- 
quency lul which is slightly detuned from the atomic 
transition frequency w. Note that for simplicity in this 
paper we identify 7r/2-pulses with Hadamard gates which 
has no effect on the estimation uncertainty. Each atom 
then undergoes a free evolution of duration t before a sec- 
ond 7r/2-pulse (using the same laser) and a measurement 
of the atomic state is performed. During the time t the 
atoms gather up a relative phase Lp = {uj — uji,)t which 
can be estimated from the measurement data. If wl and 
t are known, we therefore obtain an estimate oJest of the 
frequency uj with an uncertainty given by pol [2l| 
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which, for unbiased estimators, is simply the standard 
deviation. The uncertainty, or precision. Aw is bounded 
from below by the (quantum) Cramer- Rao bound [20l - [2^ 
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where F is the Fisher information and Fq is the quantum 
Fisher information (QFI). Expressions for F and Fq can 
be found in [Ml, [11] and Appendices [B O The Fisher in- 
formation depends on the state of the system before the 
measurement and the measurement itself while the QFI 
depends only on the state before the measurement. The 
first bound in Eq. ([2]) can be reached via maximum likeli- 
hood estimation for large ly (or T) and the second bound 
by an optimal measurement which always exists [2^ . 

If we assume that our pure input state remains pure, 
a product state iV'f™) a-s input then leads to the SQL 
precision Aw,„i„ = 1/vTtN, whereas an entangled 
(iV-particle) Greenberger-Horne-Zeilinger (GHZ) state, 
\i'?n"^) = (|0)^^ + |l)®^)/\/2, improves the precision 
to Awmin = l/VrtN, i.e. the Hcisenberg limit [l^ . 
However, under realistic conditions, the pure input state 
will degrade into a mixture due to unavoidable noise. 
For the systems considered here the dominant source of 
noise is dephasing caused by fluctuating (stray) fields 
leading to random energy shifts of the atomic levels. As 
shown in [l^, [23| , the advantage of a GHZ state deterio- 
rates in case of uncorrelated dephasing, leading to exactly 
the same optimal precision as a product state which has 
merely SQL scaling. However, this is not the situation 
commonly encountered in ion traps or atoms in optical 
lattices where particles are very closely spaced. Here, the 
particles are subject to the same fluctuations which leads 
to correlated dephasing such that the time evolution of 
the system state p is determined by 
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FIG. 2: (Color online) (a) Precision AiUmin versus duration of 
the free evolution for N = 6 atoms. The minima of the curves 
define Actf^n- (b) Best possible precision Ao;^*^ versus atom 
number A''. In both figures the dotted (red) line corresponds 
to a GHZ state, the dashed (black) line to a product state and 
the solid (blue) line to the optimal precision. In (b) we also 
show the precision corresponding to a product/GHZ state un- 
dergoing uncorrelated dephasing (dashed-dotted, green line). 
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Wi, 7 is a dephasing rate and L ~ Sz 



^j^i (^i, where cr^ is the Pauli z-operator acting on atom 
j. The fact that Eq. ([3]) describes the dominant source 
of noise in the setups considered in this paper was very 
clearly shown in a number of recent experiments [ij- 
[l3|. Equation © can be derived via a Langevin-equation 
approach by assuming that the atoms are subject to level 
shifts caused by the same fluctuating field with a sharp 
time correlation function (see Appendix |X] for details) . 
This is typically the case in experiments, e.g. in the ion 
trap experiment [l^ , where the dominant source of noise 
is caused by fluctuations of the homogeneous magnetic 
field which is required to lift Zeeman degeneracies and 
which affects all ions in an equal manner. 

Suppose we consider states which are symmetric under 
particle exchange, then we can use a Fock representation 
in which Sz = uq — ni, where rii = ajoi, and a| (a^) are 
bosonic creation (annihilation) operators of an atom in 
state i = 0, 1. We can then rewrite Eq. ^ to obtain 
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and a symmetric, pure input-state has the form Itpin) = 
J2k=o '^k\k, N — k), where \k, N ~ k) is a Fock state with 
k {N — k) atoms in state |0) (|1)). Equation (|4]) can 
be solved analytically [see Eq. (jBlSP ] which yields the 
system state immediately before the measurement which 
can be used to calculate the QFI. A GHZ state is in 
this representation formally equivalent to a NOON-state 
known from optical interferometry [3], iV'in) = (|A^jO) -I- 
|0, A^) )/ -\/2. Using this state as initial state leads, via the 
QFI, to the precision 
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The quantity Aw^j„ is found by using an optimal time 
topt = 1/'2.^N'^ for each experimental run. As can be seen. 
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^^min dependency on N and therefore there is no 

advantage using a GHZ state in the presence of cohective 
dephasing. The precision Awmm and Aa;^*„ correspond- 
ing to a product state can be calculated numerically lead- 
ing to Acj;^*„ « (^/2 -f 0.87/iV°-90)y777^ which has no 
SQL scaling and docs not even approach zero for large 
N, but is still better than the GHZ case (see Fig. (2). 

A decisive feature of collective dephasing is the exis- 
tence of decoherence free subspaces (DFSs) [l^ which 
are given by states such that ijV'DFs) = 0. However, 
since L = 5*2 in Eq. ^ a highly robust DFS state would 
be stationary and hence useless for frequency estimation. 
Ideally, one would therefore use input states which lead 
to an optimal trade-off between gain in precision and ro- 
bustness which can be found by maximizing the QFI with 
respect to all possible input states. In Appendix [B] we 
show that the maximal QFI can be attained by states 
which are symmetric under particle exchange leading to 
a considerable simplification of the optimization prob- 
lem. Results are shown in Fig. [5] As can be seen the 
best possible state leads only to a marginal improvement 
over a product state. For comparison, in Fig. [IJb) wc 
also plot the precision corresponding to iV'f™) which is 
subject to uncorrelated dephasing [13[. Evidently, cor- 
related dephasing is significantly more destructive than 
uncorrelated dephasing. 

To make use of the coherence preserving features of 
DFSs we have to alter the dynamics of the system such 
that the incoherent part of Eq. ^ is zero and the coher- 
ent part is non-zero. To this end we consider a scheme 
consisting of atoms [N even), where half of the atoms 
accumulate a phase ipi and the other half (f2- Such a 
scheme is realized in a system where, e.g., every odd atom 
has a transition frequency ui and every even atom has 
a transition frequency u)2 [see Fig. [IJb)] and our goal is 
to estimate the frequency difference 5 = uji — u)2 (see 
Appendix lA 2[) . If fiuctuating fields lead to the same 
energy shift in both transitions the incoherent part of 
the master equation ^ vanishes if a DFS state of the 
form \i!,r,) = (I0101...01) -I- |1010...10))/V2 is used. 
Via the QFI, we then obtain the bound for the precision 
A(S„i,;„ = 2/^/TiN which has Heisenberg scaling even in 
the presence of correlated dephasing. We should note 
here that in general we can use arbitrary orderings of the 
atoms in Fig.[TJb). The input state then takes the form 

1 ^ 

\'4'tn) = — =(|il,l2, • ■ • , iw) + J]^ 0'^|jl,«2, ■ • ■ ,iAr)), (6) 

where ij = 0, 1 and X^jLo *i ~ (i-C- ij ~ occurs as 
many times as ij ~ 1) and atoms with ij — (ij ~ 1) ac- 
cumulate the phase iy9i {(p2)- The described dynamics can 
be obtained, e.g., by choosing the two transitions to be 
within the same Zeeman manifold such that the difference 
of the magnetic quantum numbers of each transition is 
equal, and thus a (weak) fluctuating magnetic field leads 
to the same energy shifts. This was demonstrated in a 
recent experiment with two ions in a linear Paul trap re- 



vealing a significant increase in the coherence time jl5l |. 
The same ideas were furthered by an experimental study 
of non-perfect input states [25| . In both experiments, an 
additional electric quadrupole field was used to obtain 
LUi L^2 and from the measured frequency difference the 
electric quadrupole moment was determined. 

The above experiments also offers an alternative view 
of the fact that Heisenberg scaling can be obtained in 
this setup. In Ref. [l^ a 'designer atom' was constructed 
consisting of two physical atoms with two internal, logical 
states |0)l = |01) and = |10) which are decoherence 
free. A state of the form (|0101 ... 01) -h |1010 . . . 10))/ 
is then equivalent to a decoherence free GHZ state 
of 71 = N 12 designer atoms, |V'm)L = (|00...0)l -|- 
|11 . . . l)i)/-\/2. The states |0)l, |1)l accumulate a rel- 
ative phase and so it is straightforward that |V'm)L 
leads to a sensitivity lS.8rain ~ l/'s/Ttn. 

We can also conceive a situation where the fluctuat- 
ing field shifts the transition frequency wi of half of the 
atoms and the transition frequency uj2 of the other half 
by the same magnitude but opposite sign. For the setup 
shown in Fig. IHb), this means that we have to replace 

the noise operator in Eq. ([3|) by L ~ Tlij=ii~^)'' '^z and a 
GHZ state would be decoherence free. Wc can utilize this 
for quantum enhanced precision measurements by per- 
forming a Ramsey-type experiment but now with up to 
two lasers of frequency loli and ujl2 such that half of the 
atoms accumulate a relative phase ipi = [bJi—bJi^ijt = 5it 
and the other half ip2 = ("^2 — ^L2)t = &2t- The Hamil- 
tonian can then be written as (see Appendix lA 3|) 

H=-^{5i+52)S,. + -^{52-5i)L, (7) 

where the second term vanishes if applied to a GHZ 
state. If the laser frequencies are known, the quan- 
tity which can be estimated with this setup is therefore 
given by = {lji + W2)/2 and the corresponding preci- 
sion, which can be calculated via the QFI, is given by 
Afimin = l/y/rtN which has Heisenberg scaling. The 
above discussion can again be generalized to an arbitrary 
ordering of the atoms as long as we use a GHZ state as in- 
put. The setup can be realized, e.g., by using transitions 
|m) -f^ |m) with frequency wi and | — m) -H- | — m) with 
frequency W2, i-e. the two ground states (with magnetic 
quantum numbers ±m) and the two excited states (with 
magnetic quantum numbers ±m) are in the same Zeeman 
manifold, respectively, such that (fluctuating) magnetic 
fields cause first order Zeeman shifts of the same magni- 
tude but opposite sign [2^ [2^ . Note that the quantity 
to be estimated, f2, is magnetic field independent (in first 
order), i.e. the situation is similar to a clock transition, 
and might therefore serve as a frequency standard. 
Moreover, H can easily be chosen to be in the optical 
domain which is desirable for atomic clocks [ll| . 

We also note that, similar to the case of estimating S, 
we can introduce decoherence free logical states |0)l = 
|00) and |1)l = |ll) which accumulate a relative phase 
2^lt. A GHZ state of N atoms is then simply a GHZ 
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state of n = N/2 logical states and the precision is given 
by Ar2„ji„ = l/y/Ti2n, the factor of 2 arising from the 
factor of 2 in the relative phase of |0)l and 

An optimal measurement for the two schemes discussed 
above, i.e. a measurement for which F = Fq [sec Eq. ([2])], 
is given by a 7r/2-pulse and a measurement of each atom 
in the {|0), |l)}-basis [see Fig. [Hb)]. However, in prac- 
tice there will be imperfections both in the preparation of 
the input state and the measurement. A faulty measure- 
ment of an atom can be modeled using the measurement 
operators 

n. = i(i + 77A/)|j>(*l + ^(i-?7MKN>(«k,, (8) 

where i = 0, 1, and i]m is the likeliness that we get the 
correct measurement result. Also, we assume that the 
actual input state is of the form 

p,;„ = ^(iV)|V',;„)(V'™| + - ^W]t, (9) 

i.e. we prepare the ideal input state |V'm) with fidelity 

/ = {4'^n\p^n\4'^n) = ^N) + ' ^W]. (10) 

and hence / « ^{N) for > 1. The error model (O is 
a worst case scenario since the identity matrix does not 
yield any phase information. Analogously, we assume 
that a 7r/2-pulse is given by the operation 
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i.e. 'qn characterizes the probability to perform a perfect 
7r/2-pulse leading to an ideal state pid (see Appendix [C|). 
We can then calculate the Fisher information for both 
schemes and therefore the Cramer-Rao bounds. 
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(12) 



where we assumed that tpi ± (^2 = t^/N ('-|-' for Afi,„i„; 
'-' for A^min), which can always be achieved by a feed- 
back setup which appropriately adjusts, e.g., the elec- 
tric quadrupole field, laser frequencies or/and the evo- 
lution time t. We note that the state which was pre- 
pared in [2^ leads to the same result for £,{N) = 1/2 and 
N = 2. The term CiN)r]^r]^i in Eq. ^ might, at first 
glance, lead to the conclusion that faulty state prepa- 
ration and detection annihilates the advantage gained 
by using a DFS. To show that this is not the case 
we compare the precisions (fT2)) to those obtained us- 
ing conventional Ramsey spectroscopy with a product- 
state I'f/'f™) as input. In this case we would use A^/2 
atoms to estimate wi and the others to estimate W2- 
For a fair comparison we assume uncorrelated dephas- 
ing which can in principle always be achieved by plac- 
ing the atoms in different traps. Assuming that \ipf^°) 
is created by N 7r/2-pulses, the corresponding precision 
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FIG. 3: (Color online) (a) Minimum fidelity of the input state 
versus number of atoms N [cf. Eq. (|13[)]. (b) Maximum like- 
lihood estimation uncertainty Ail versus total measurement 
time T for iV = 20 and ^(iV) = 0.6. The lower, solid (blue) 
line corresponds to a GHZ state and the upper, solid (red) 
line to a product state. The dashed (black) line is given by 
Eq. ini). In both figures we set r]H = 0.98, rjM = 0.99, ■yt = 3. 



then reads AuJi^min = Aw2, mm = ^/■^cy /NT /rjjjrjM , and 
the precisions of the quantities to be estimated is given 
by Afimin = \l\5„yin ~ Aa;i,,„i„/V2 which have to be 
compared to Eq. ([T2|) . In both cases this leads to the 
constraint 
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2 N~ 



(13) 



As can be seen the bound for the state fidelity £,{N]_ is 
surprisingly low. In the experiment described 



i.e. whenever the above inequality is fulfilled the DFS 
schemes beat conventional Ramsey spectroscopy. An ex- 
ample is shown in Fig.|3l^a). With current ion trap exper- 
iments gate and readout fidelities in excess of r/n — 0.98 
and iiM = 0.99 have been achieved Furthermore, we 
assumed that using a DFS scheme leads to a coherence 
time which is 3 times longer than the coherence time of 
a single atom. This is a rather conservative estimate 
which has already been exceeded in experiments [l3 |. 

in [3 a 

50.8% fidelity for a A^ = 14 GHZ state was achieved. For 
our scheme it would be required to manipulate this GHZ 
state by transfcring half of the atoms into a different in- 
ternal state. Naturally this would be done by addressing, 
for example, the second N/2 neighbouring atoms by an 
appropriate sequence of laser pulses, i.e., crucially, it is 
not required to address atoms individually. This would 
of course decrease ^{N) but even if it reduces it to, say, 
20% (which is a very conservative assumption) we still 
beat conventional Ramsey spectroscopy. 

The bounds (jl2p can be reached using maximum like- 
lihood estimation in the limit of large u (or T). In prac- 
tice it is certainly highly relevant how large v has to be 
such that the actual estimation uncertainty is close to 
the bound. Suppose we perform v experimental runs and 



obtain the results 



ni. 



where n,- is the number of 



times the state |0) is measured in each run, and the total 
number of even rij is Ve- It turns out that the maxi- 
mum likelihood estimators flest and (Jest depend only on 
(sec Appendix [D|. Using the probability distribution 
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for Ve and Eq. ([T]) we can then calculate, e.g. the esti- 
mation uncertainty Af2 for finite v. A result is shown in 
Fig. El^b) depending on the total time T = vt oi the ex- 
periment (lower solid line) for (pi + if2 = t^/N. As can be 
seen the estimation uncertainty quickly approaches the 
lower bound (dashed line). We also show the estimation 
uncertainty and the lower bound for conventional Ram- 
sey spectroscopy with ji/'f™) (upper solid line; the two 
quantities are indistinguishable on the scale of the fig- 
ure). Evidently, even for small T, the DFS scheme easily 
outperforms conventional Ramsey spectroscopy. We note 
that the corresponding plots for estimating 6 would be 
identical to the ones shown but larger by a factor of two. 

To conclude, we have shown that correlated dephasing 
significantly diminishes the precision of frequency estima- 
tion with standard Ramsey interferometry. On the other 
hand, it allows for the existence of DFSs which we used to 
construct and analyze generalized Ramsey setups which 
beat the SQL even in the presence of faulty detection and 
significantly imperfect state preparation. The proposed 
schemes for quantum enhanced frequency estimation are 
therefore feasible with current experimental technology 
and can lead to improved spectroscopic methods with a 
variety of important applications in metrology. 
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Appendix A: Noise model 

In this appendix we give a detailed description of the 
noise model and the derivation of Eq. ([3]). An alter- 
native derivation is given by coupling the atoms to a 
bosonic bath, similar to the methods described in (28l . [29j . 
However, the following derivation, which is based on a 
Langevin-equation approach, is physically more intuitive 
for the systems considered in this paper, i.e. atoms which 
are subject to fluctuating classical fields. 

Consider N two-level atoms (or ions) with internal 
states {|0)j, and transition frequencies ujj which are 
subject to a time-dependent, fluctuating field leading to 
random energy shifts of the transitions. The Hamiltonian 
can then be written as 



N 



H = 



Ho + B{t)L, 



N 



(Al) 



where (the real numbers) B(t) and Ej characterize the 
field strength and how it affects atom j, and cr^ = 
|0)j(0| — is the Pauli z-operator acting on atom 

j. Writing i?(i) = /2^{t) the Schrodinger equation 



then takes the form 



-iHom-i^l^Lmit)- 



(A2) 



The random fluctuations of the field are captured in ^ [t) 
and Eg. (jA2| is therefore an example of a Langevin equa- 
tion [30| . Assuming that ^(t) has zero mean and very 
rapidly decaying time correlations, we can write 



m - 0, 



mm = 5{t~t'), 



(A3) 
(A4) 



where the overbar denotes the mean value, i.e. we make 
the idealization that £,{t) is white noise. Equation (|A2p 
can then be written as a stochastic differential equation 



(5), (A5) 



where dW = ^{t)dt is a Wiener increment [301 . The 
(S) indicates that we have to interpret this equation in 
the Stratonovich sense. The reason for this is given by 
the fact that Eq. (jA4[) is an idealization. In reality this 
correlation function will have a finite width (but which is 
small compared to any other relevant time scale) . In such 
a case the stochastic differential equation (|A5p has to be 
interpreted in the Stratonovich sense. A detailed discus- 
sion of this point can be found in Chapter 6.5 in [30j . 
Transforming Eq. (|A5|) into Ito form, using the standard 
rules [301, leads to 



m = {-iHo - \^)dt 



Ul^L\i;)dW (/), 



(A6) 

where the (/) indicates that this equation is to be in- 
terpreted in the Ito sense. Using Ito calculus, the time 
evolution for the density operator g= \ip){'ijj\ is then de- 
rived to be 

dg = g{t + dt) — g{t) 

= m+dt)){^{t+dt)\-\m){m\ 
= {m)) + \dmm)\ + m)~\m)m)\ 
= \di,){^\ + \^){d^\ + \di:){d^ 

= (-ii?o - Ji^) gdt + g (ii/o - Ji^) dt 



-ul'^LgdW + iJ'^gLdW +^LgLdt 



For the averaged density operator 

pit) - m 

we therefore get 



p = -i[Ho,p] 



^ LpL ^r2^ Kr2 



2 V 



-L^p- -pU 
2 ' 2^ 



{!)■ 
(A7) 

(A8) 
(A9) 



In this paper we consider three different scenarios corre- 
sponding to special cases of the above equation, some of 
them make use of dccohcrence free subspaccs (DFSs): 
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1. Conventional Ramsey Spectroscopy 

In conventional Ramsey spectroscopy we consider the 
same transition in each atom, i.e. lo = ujj, j = 1 . . . N 
and hence all transitions will be affected in the same way 
by the fluctuating field, i.e. £j = 1, j = 1 . . . N. It follows 
that we have (in a rotating frame with respect to a laser 
frequency wl) 



N , _ s 

L = S.^Y.<^ Ho=^^^^S., (AlO) 
and thus Eq. (|A9|) is equal to Eq. jS]). 



frame with respect to laser frequencies uj^i and iOL2) 

L = -T.< + 11< (A13) 



2 2^'^^+ 2 



'jeA 



= '-^S. + ^L, (A14) 
where Si = oji — utLi- An iV-particle GHZ state 

|^™) = i=(|00...0) + 111...1)) (A15) 

has the property L\ipin) ~ 0. This setup can be used to 
estimate 51 = [coi +a;2)/2. 



2. DFS spectroscopy for estimating difference of 
two frequencies. 



Here we assume that half of the atoms, represented by 
a set A, have transition frequency uii and the other half, 
represented by a set B, have transition frequency i02 and 
the fluctuating fleld leads to the same energy shift in both 
transitions, i.e. = 1, j = 1 . . . A^. Hence we obtain 



L 
Ho 



5., 



UJi — UJ2 



Ull + IjJ2 



L.(All) 



If we use a state of the form 
1 



T2' 



N 

ii,i2, . . . ,ijv) + ]^cr;J.|ii,i2,---,iAf)), 



(A12) 

where ij = if j G A, ij = 1 if j g B and = \0)j{l\ + 
|l)j(0|, we have L\ipin) = and the above setup can be 
used for the estimation of (5 = cji — uj2- 



3. DFS spectroscopy for estimating mean of two 
frequencies. 



Finally, we assume that half of the atoms, represented 
by a set A, have transition frequency uji and the other 
half, represented by a set B, have transition frequency 
UJ2 and the fluctuating field leads to energy shifts of the 
two transitions of the same magnitude but opposite sign. 
More exactly, ivj = oji, ej = — 1 if j € A and ojj = 
^2, £j = 1 if J G -B. In this case we obtain (in a rotating 



Appendix B: Quantum Fisher information 

Consider a system state p which depends on a param- 
eter a which is to be estimated. Defining p' = j^p, the 
quantum Fisher information (QFI) is given by 



Fq = Tr{p'Cp{p')}, 



(Bl) 



where Cp{p') is the "symmetric logarithmic derivative" 
(SLD) of p H^^ii, mf. Writing the state in diagonal 
form, p ~ X]j iV'j) (^j |i the SLD is given by 



E 



and therefore 



Pj + Pk 



{^,\p'\^PkM,){yJk\ (B2) 



E 



P] + Pk 



(B3) 



Assume that p is the solution of Eq. (|A9p and [L, Hq] = 
and that only Hq depends on the parameter a such 
that [Hi), Ho] = 0, where i/p = j^Ho- We then obtain 
p = e~'^"*pe'-^"*, where p is the solution of Eq. (jA9P with 
Hq = 0, and therefore 



^^Q = 2t^E 



J,k 



jPj -Pkf 
P] + Pk 



\mHM? 



(B4) 



Note that if p = |'0)(V'l is a pure state, the above reduces 

toFQ = Aem{H'QfH)-{m'oW)- 

In conventional Ramsey spectroscopy (see Ap- 
pendix lA II) a GHZ state of the form (jAlSp would evolve 
into 



Pit) 



1 r 

2 



|0...0)(0...0| + |1...1)(1...1| 



(e-"5^*|0 . . . 0) (1 . . . 1| + e'*^*|l ... 1) (0 ... 0|) 



(B5) 
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where 5 = lo — ujl. As can be seen, due to corre- 
lated dephasing, the above state decoheres on a timescale 
which is shorter than the decoherence timescale 
in the presence of uncorrelated dephasing (given by 
l/'jN). This behavior was therefore dubbed "superde- 
coherence" [3, [11] . The parameter to be estimated is 
the transition frequency a = w. The corresponding QFI 
is obtained by diagonalizing the state (|B5|) and using 
Eq. (HH leading to Fq = t'^N'^c''^^'^^'^ . The correspond- 
ing precision 



AuJn 



1 



(B6) 



opt 



is optimal for a time topt = 1/2'yN'^ leading to A; 
y^2ej/T which, as a consequence of superdecoherence, 
has no N dependency. GHZ states are therefore not par- 
ticularly useful for conventional Ramsey spectroscopy. 

To find the best possible precision in conventional 
Ramsey interferometry we can restrict ourselves to in- 
put states which are symmetric under exchange of parti- 
cles. To prove this we use a method inspired by Ref. [s^ . 
Consider a unitary operation U which maps an arbitrary 
state onto a symmetric state, i.e. 



N 



U\i>)=J2ck\k,N-k). 



(B7) 



k=0 



The state \k,N k) is the completely symmetrized state 
with k atoms in state |0) and N — k atoms in state |1) 
and is defined by 



\k,N -k) = 




11,^2, ... ,i7v), (B8) 



where ij = 0, 1 and k (N — k) is the number of zeros 
(ones) in 12, . . . , jtv)- Furthermore, the sum is over 
all permutations P of particles which lead to different 
terms in the sum. Before we proceed we will show (i) 
that such an U always exists and (m) that [J7, Sz] = 0. 

Proof of (i): An arbitrary state of the system can be 
written as 



N Mk 



(B9) 



where |fc;/^) = |ii,i2, 



fc=0 /j=l 



0, 1 and k is the num- 



ber of times ij is zero. For each k there are Mk — v fe / 
such states which we enumerate using /i. On the sub- 
space defined by a fixed k we define a basis {|0j(fc)) \i ~ 
such that 101 (fc)) = VmI^;M> with 

= \JY.^=i l^fc.^P and the remaining {|(/'j(fc)) | j > 1} 
are chosen such that we obtain an orthonormal basis. 
Equation (lB7| is fulfilled if U\(j)i{k)) = jf^\k,N - k) 
for all k. Decomposing U into a block diagonal form 
U = ©^^g U{k), i.e. U{k) is the Mk x Mk block acting on 



the subspace fc, the above can be achieved by requiring 
that the first column of the matrix U{k) has elements 
U,,{k) = {4>,{k)\U{k)\Mk)) = f-^{^,{k)\k.N - k), 
j = l,...,Mfc. The remaining columns can be chosen 
such that all columns are mutually orthonormal. If A/fc 
is zero for some k we can choose U (fc) to be the identity. 
Therefore U{k) and hence U can be chosen to be unitary. 
Proof of (ii): We have 

JV 



SzU\i^)^SzY.Ck\KN ~k) 

k=0 

N 

= Y^{2k- N)ck\k,N -k) 



k=0 

N Mk 

= Y,i^k - N)U ^bk.Jk; ij) = USz\^) 

k=0 M=l 

(BIO) 

and hence [U, Sz] = 0. 

We can now use (i) and {ii) to show that to find the 
optimal precision we can restrict ourselves to the sym- 
metric subspace. Consider an arbitrary pure input state 
\4^in)- The time evolution of this state in conventional 
Ramsey interferometry as defined in Appendix lA II and 
Eq. ([3]) is given by 



Pit) =e-'^^^*e-i^-'* 
(7t/2)- 



E 

m=0 



■5ri0m)(^m|5re--^^*e'- 



(Bll) 



If we take the symmetric state U\ipin) StS input state 
instead of iV'in) the state of the system at time t has 
the form p*(<) = Up{t)W due to {ii). Diagonalizing 
p{t) using an orthonormal basis we can write p{t) = 
Y.kPk\'^k){^k \ and therefore p''{t) = Y.kPk\i'k){4'k\ with 
1^^.) = U\ipk)- Due to {ii) the QFIs of p" (t) and p{t) are 
therefore equal, 



71 P^+P'^ 



Pj + Pk 



mszmi^ ^ fqIp]. 



(B12) 



Assuming that {tpin) is an optimal input state which max- 
imizes the QFI then U\^in) is optimal as well and there- 
fore we can restrict our search to the symmetric subspace 
which concludes the proof. 

The Fock states defined by Eq. (|B8|) represent states 
with k atoms in state |0) and N — k atoms in state 
Also the operator Sz can be written in a Fock representa- 
tion given by Sz ~ uq — ni, where rii ~ aj^a; and Oi (aj) 
arc bosonic annihilation (creation) operators for modes 
i = 0,1. Since the total particle number no + ni = N is 
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conserved we can set 5*2 = 2no — N and therefore Eq. ([3]) 
transforms into 

p = -i(5[?io, p] + 27 (^nopno - hnlp - ^P?^o^ • (B13) 

Furthermore, every symmetric, pure input state can be 
written in the form 



N 



k=0 



The solution of Eq. (|B13[) is given by 



Pit) 



E 

m=0 



-n^p{Q)n^e 



(B14) 



N 



afeare-T*('=-'' e-'^'^^'-'^l/fc, - k){l,N ~ II 

(B15) 



where we set p{0) = \ipin){^in\- In order to find the 
input state which leads to the best possible precision for 
estimating uj we performed a numerical optimization in 
the bosonic picture, where the time evolution is given by 
Eq. ()B13|) using methods described in [l^l and a result is 
shown in Fig. [2jb). 
A product state 



IC°> 



^^(|0> + |1)) 



(B16) 



is symmetric under particle exchange, and in the Fock 
representation it takes the form of a 'coherent state' 



1 ^ 

^ ^ fc=0 




\k,N -k) 



(4 + 4)^10,0). 



(B17) 



Using this as input state, the density matrix (jB15[) can 
be numerically diagonalized, and via Eq. (|B4I) we cal- 
culate the QFI. Like for a GHZ state, Auj„iin can be 
minimized for a time top* and the corresponding preci- 
sion is obtained to be Aa;°^*^ « {^2 + 0.87 /N^■^°)^/r/T■ 
This shows that, also in the case of product states, corre- 
lated dephasing is more detrimental than uncorrelated 
dephasing (in which case we would obtain Au)'^*^ = 

^^ejN ^ ^ jT [3]). Furthermore, the best possible pre- 
cision is only marginally better than the precision ob- 
tained by using a product state [see Fig. [2l^b)] showing 
that conventional Ramsey spectroscopy is merely of lim- 
ited use for frequency estimation in the presence of col- 
lective dephasing. 



Appendix C: Fisher information 

The QFI provides the optimal precision for estimating 
a parameter. It depends only of the system state before 



the measurement and not on the measurement itself. In 
order to examine the effects of particular measurements 
on the estimation precision we therefore have to consider 
the Fisher information (FI). The FI is given by 



^ p{k\oL) \da 



(CI) 



where p{k\a) is the probability to obtain a measurement 
outcome k given that the value of the parameter to be 
estimated is a, 



p{k\a) = Tr{Ukpia)}. 



(C2) 



Here, the operators Ilk form a positive operator valued 
measure (POVM) describing the measurement. If for a 
particular POVM the FI is equal to the QFI the measure- 
ment is said to be optimal, i.e. it saturates the quantum 
Cramer- Rao bound [see Eq. ([2])]. 

Both for the estimation of a = 5 = cji — a;2 and 
a = U = {uji + 012)72, i.e. the two schemes described 
in Appendices lA 21 and lA 31 the optimal measurement is 
given by a measurement of all atoms in the cr^; -basis which 
in practice is done by a Hadamard gate and a measure- 
ment in the {|0), |l)}-basis. Note that we use Hadamard 
gates for simplicity. In practice these can be replaced by 
7r/2-pulses which has no effect on the FI. In an actual 
experiment both Hadamard gate and measurement will 
have imperfections. To model these we assume that the 
Hadamard operation on one atom is given by 

£Hip)^VHHpH+^{l-7^H)t, (C3) 

where iJ is a perfect Hadamard gate and tjh characterizes 
the probability to have a perfect gate. It corresponds to 
the gate fidelity fn, as defined e.g. in [13], via Jh = 
■\/(l + i1h)/2. The POVM for the measurement of one 
atom is given by 

no = ^|o)(o| + i^|i)(i|. 



Hi 



2 

1 + Vm 



2 

1 - Vm 



(C4) 



i.e. r]M quantifies the probability that we have a perfect 
measurement. The above can be combined into a new 
POVM which describes a faulty measurement in the ax- 
basis, 

n. = i±f^|±)(±| + i^f^|T)(T|, (C5) 

where |±) = (|0) ± |l))/\/2 arc eigenstates of ctx- Imper- 
fect state preparation can be modeled by 

P^n = aN)\An){i^^n\ + ^ [1 " ^(A^)]!, (C6) 

where {ipin) is the ideal, pure input state. For the es- 
timation oi a = 5 = uJi — UJ2 (see Appendix lA 2[) the 
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state \tjjin) is given by Eq. (|A12p and for the estimation 
of a = = (wi +a;2)/2 (see Appendix lA 3p the state 
lipin) is given by Eq. (jAlSj) . The state pin evolves then 
into the state p{a), the state before the measurement, 
according to the dynamics given by the Hamihonians in 
Appendix IA2I and Appendix \A3\ respectively. 

A particular outcome fc of a measurement on all N 
atoms is given by a sequence {11,12,13, ■■■ ,iN} where 



±, i.e. if the jth atom is found in state 



(h)) 



we have ij = + {ij = —). Note that in practice a |+) 
(|— )) outcome corresponds to finding the atoms in state 
|0) (|1)) due to the Hadamard gate. The probability for 
a particular outcome is then calculated to be 

p(fc|a) = Tr{Tli^Ili^Ili^ . . .Ili^p{a)} 



1 



(1 + i~miN)v^iv^i cos[N^ia)]) 



2N 

Quia), 



where 



(p{a) 




ujL2)/2]t ; if a = 17 
; if Q = (5, 



(C7) 



(C8) 



and n is the number of times '+' is contained in the 
sequence {ii, 12, 13, ... , in}- From this we obtain the PI 



F 



N 

E 



1 



n ) qn (a) \ da 



■J-Qnia) 



c(a)jVte(A^)?ygO sin^[iVy(a)] 
l-(CW«)'cos2[]V^(a)] 



iC9) 

where c(a = 6) = 1/2 and c(a = = 1. The FI is 
maximized for ip{a) ~ tt/2N leading to the Cramer-Rao 
bound 



(CIO) 



2 ^tN£,{N)^^TiZ' 

Setting V — T/t wc therefore obtain Eq. ([T2l) . 

Equation (jC10| has to be compared to the precision 
corresponding to 'classical' Ramsey spectroscopy using 
the product state (|B16[) as input, i.e. a scheme which 
does not rely on non-classical correlations between the 
atoms. Wc will assume in the following that in this case 
the system is subject to uncorrelated dephasing since the 
atoms can in principle always be put in separate setups. 
To estimate VI or 5 with this method we use A^/2 atoms, 
represented by a set A, to estimate wi and the remaining 
A^/2 atoms, represented by a set i?, to estimate UJ2- Since 
the system state is a product state these two estimations 
are completely independent. 

Preparation of the state (jB16p is achieved by 
Hadamard gates e.g. on the state |00 . . . 0), i.e. the state 
of an atom j before the measurement is given by 

P.(i)-U|0)(0| + |l)(l| 



-7* 



'|0)(1| 



|1)(0|) , (Cll) 



where ej = cji — oj^i = ea for j £ A and ej = 0^2 — i^L2 = 
Eb for j ^ B. The atoms are measured in the a^j-basis 
described by the POVM ()C5p . and therefore we obtain, 
e.g. for the atoms in group A 



P 



(fc|c.l) = P+(c.l)"P_(c.l)^/2- 



(C12) 



where 



P±(c^i) = - (1 ±ry2^77Me-^*cos(eAt)) , (C13) 

and n is the number of measurement outcomes of the 
atoms in group A. With the help of this we obtain 

•^l^iT^l^Me-r^fcos^EAt) 

which is maximal for ea ~ """/St, and for t = topt — I/27 
we have 



ije 1 



(C15) 



The corresponding expression for Au)2.min is obviously 
the same and therefore we have 



An, 



A(5mm AoJi^. 



V2 



(C16) 



The above expression has to be compared with Eq. (|C10[) 
leading to Eq. ([1] 



Appendix D: Maximum likelihood estimation 

To construct the Maximum likelihood estimators fiest 
and Sest corresponding to the DPS schemes described in 
Appendices lA 21 and lA 3[ we consider a sequence of i' ex- 
perimental runs with results ni , . . . , , where Uj is the 
number of times we obtain the result '+' in the jth rep- 
etition of the experiment. Using (a) from Eq. (jC7p . 
the likelihood function for such an outcome is given by 



£,{a\ni,...,n^) = Y[[^ 



X {l^aN)v^ivZcos[N^{a)]Y- 



1 -TT rN 



n 



(Dl) 



where <p(a) is given by Eq. (jCSp . and Ve is the number 
of even rij . Maximizing Eq. (jPip with respect to a leads 
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to the estimators 



Wt ai-CCOS ( ^^{N),jN^ 



Nt 





< I^e < l^i 
■,Ve<V- 



(D2) 



where v± ~ v{l±S^{N)r]^ri^.j) /2. The probabihty distri- 
bution for v„ is calculated to be 



(l + eW^Mcos[iV^(a)] 
(l"CW'7Mcos[iV^(a)] 



(D3) 



where (p{a) is again given by Eq. (jCSp . Using p(!^e) we 
can numerically calculate the first and second moments 
of ilesi and 5est leading to the precision of estimating, 
e.g., for finite v = T jt^ 



An 



^d{n,,t)/dn\ 




(D4) 



which is shown in Fig.[3jb). 

For the sake of completeness we also give the maxi- 
mum likelihood estimator if the input state is a product 



state (jBlSp and undergoes uncorrelated dephasing as dis- 
cussed in Sec. [C] Like before, the goal is the estimation 
of u>i and 0J2 using N/2 atoms for each, from which we 
can estimate fl and S. Since neither the state nor the 
noise nor the measurement are correlated we can treat 
the problem of, e.g., estimating tui as if there is only 
one atom which is measured i> = Ni>/2 times with pos- 
sible outcomes ij = ±. Denoting the total number of 
'+'-outcomes n, we obtain 



(D5) 



where P± {uji ) is given by Eq. (|C13p . Maximising £ leads 
to 



Widest ^ '^Ll 



.0 



\NvrijjriMC ^* ^ ' — — + 

■n<u_ (D6) 



where v± ~ Nv{l ± rjjjrjMe '''*)/4. The probability dis- 
tribution for n is 



(D7) 



The corresponding expressions for estimating W2 are of 
course identical. The estimators for Q, and 5 are then 
simply given by Vt^st = (1^1, est + i^2,est)/2 and Jest = 
wi,est —^2, est- Usiug Eqs. (|D6p and (|D7p we can numeri- 
cally calculate Aft and A6 for finite ly = T/t. The former 
is shown in Fig.[3][b). 
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